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Abstract 

In the framework of finite temperature conformal scalar field theory on de Sitter 
space-time the linearized Einstein equations for the renormalized stress tensor are 
exactly solved. In this theory quantum field fluctuations are concentrated near two 
spheres of the de Sitter radius, propagating as light wave fronts. Related cosmo- 
logical aspects are shortly discussed. The analysis, performed for flat expanding 
universe, shows exponential damping of the back-reaction effects far from these 
spherical objects. The obtained solutions for the semiclassical Einstein equations in 
de Sitter background can be straightforwardly extended also to the anti-de Sitter 
geometry. 
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1 Introduction 



In the modern picture of the early universe the exponential expansion period is known 
to be successful in solving a number of important problems occurring in the standard 
cosmology. During this inflationary phase, described by de Sitter geometry, quantum 
field effects played an essential role for subsequent formation of the present universe. 
In this framework, the knowledge of the corresponding quantum state, originating from 
Planck epoch, is important to reconstruct the phase transitions pattern of the inflation. 
One of the possibilities generally used JI|], is to consider this quantum state as a de 
Sitter invariant vacuum experienced by a freely moving observer as a thermal bath at the 
Hawking temperature fi^ 1 0. 

A wider class of quantum states, associated in static de Sitter coordinates to a thermal 
equilibrium at arbitrary temperature /3 _1 , has been investigated in our previous work ||. 
In the present letter we extend this analysis studying the back-reaction on de Sitter 
geometry in these quantum states. 

For simplicity we restrict our consideration to massless conformally invariant scalar 
theory where the finite temperature Green function, globally defined on the whole space- 
time, can be exactly found. However, the corresponding renormalized stress tensor for 
arbitrary temperatures includes a traceless part singular on the horizon surface of two 
antipodal static coordinate systems, where thermal equilibrium is introduced. This ex- 
tremely large energy density implies that, very close to the horizon, physical quantities 
can be computed only in the framework of quantum gravity, theory which would provide 
a natural short distance cut-off []. Thus, the singular part of the stress tensor can be 
interpreted as matter distribution sharply concentrated over two spherical shells of radius 
Tq (horizon size) that propagate through the de Sitter hyperboloid as light wave fronts. So 
far as, apart from the rest space-time, the quantum field inside the domains bounded by 
the shells is in thermal equilibrium, we will refer to these spherical 3-d regions as bubbles 

0, that can be characterized by temperature (3 1 , internal energy, and entropy. 

1 Note that recently much attention has been paid to similar divergences appearing in statistical 
mechanical computations of black hole entropy when temperature different from the Hawking one 

has to be introduced to obtain the derivative of the partition function. Possible role of superstring theory, 
providing such fundamental cut-off, is discussed in 0. 

2 It is worth reminding that these bubbles are different from the topological defects nucleated during 
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This letter is organized as follows: in section 2 we compute the renormalized stress 
tensor at finite temperature (3 ^ /3h, and discuss its properties inside and outside the 
bubbles. In section 3 this tensor is used to evaluate, in the linear approximation, the back- 
reaction effects. The model of the flat expanding de Sitter universe is then considered in 
section 4 to study the cosmological implications of the chosen quantum state. Finally, we 
give our conclusions and remarks. 

2 The quantum state and stress tensor 

Let us determine the average value of the energy momentum tensor for conformally invari- 
ant scalar field theory < T£ >p= T^{f3) in static de Sitter space, where j3 is the inverse 
temperature of the system. This tensor can be decomposed as 

TZ(P)=f»(f3) + ±6;T , (2.1) 

where T"(/3) is the traceless part, and T = T^(j3) is the local conformal anomaly that 
does not depend on the quantum state of the system || . In particular, for de Sitter space 

T = (2.2) 

(r stands for the de Sitter radius), and due to the conservation law T" = 0, the traceless 
part has only one independent component. Consequently, to find all its components it is 
sufficient to calculate only one of them. 

It is worth reminding that static de Sitter metric reads 

ds 2 = g^dx^dx" = (1 - r 2 /r 2 )dt 2 - (1 - r 2 /rl)~ l dr 2 - r 2 (d6 2 + sin 2 9d V 2 ) , (2.3) 

where x M = (t,r,9,ip). In this space the temporal component of the stress tensor TJf(/3) 
can be obtained observing that [[!(] 



<f/> /3 -<f/> /3=00 =^^ . (2.4) 

■jU 9tt 



The last equation can be also rewritten in the following form 



< f! ><=< f! >* 4 f 1 - 1) S • < 2 - 5 ' 



phase transitions 



where (3h = 27rr denotes the inverse Hawking temperature. Hence, due to the fact that 
for (3 = (3h the stress tensor is completely anomalous (T"(Ph) = 0), we have 

i?<0 f 1 - Jf) £ • (2.6) 

This enables to obtain the stress tensor in the form 

< ^ > ^s( 1 -|)^ !O 4 1 '-?-5'4) + 960W^ ■ < 27 > 
that coincides at zero temperature with the known result for the static conformal vacuum 
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As one can see from ( [2.7| ) the traceless part dominates near the horizon and corresponds 
for {3 < (3h to the energy momentum tensor of a gas of massless scalar particles [|IJ| . In 
this case, the energy density has the usual Planck form and results to be 



Tm = s/" 4 ' (2 ' 8) 

where = /3 _1 (1 — Z^ 4 //^) 1 / 4 ^ 1 ^ 2 plays the role of a local, redshifted temperature. 
The similar structure for the energy density in finite temperature Rindler space is given 
TBI. 



in 



Let us point out that ( |2.7|) can be obtained by the standard procedure from the finite 
temperature Green function given by JTTJ] 

PK ' 1 2/fe/3sinha 1 cosh(ai/3 H //3) - cosh[(t - f)/3 H /r Q /3] ' 1 ' 1 

(coshai = [(r-Q — T 2 )(rQ — r /2 )]~ 1 / 2 {r 2 — rr'fcos 9 cos 9' + sin 9 sm9' cos(y? — v 5 ')]}); with the 
subtraction of its value corresponding to the Hawking temperature. Although expression 
( |2.9| ) is defined only in the region bounded by the horizon surface, it can be extended 
to the total de Sitter space-time. In terms of conformal diagram, figure 1, where static 
coordinate system ( |2.3| ) maps into the region (a), it is easy to prove that G^x, x') can be 
smoothly continued from the region (a) into the casually connected domains (c) or (d), 
by moving one by one its arguments through the horizon surface. Then extending this 
procedure from (c) or (d) to (b) we obtain a globally defined two-point Green function, 
smooth everywhere for not coinciding arguments. Such procedure completely fixes in all 
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the de Sitter space-time the quantum state, whose properties can be derived from the 
structure of the associated global energy momentum tensor. This global stress tensor as 
opposite to the Green function is singular on the horizon surface, which means that close 
to the horizons our semiclassical approach is not applicable and the complete theory of 
quantum gravity should be used for computing the quantum effects. 

As we have already mentioned one can interpret the energy singularity as a matter 
distribution over the surface of two spherical bubbles having the horizon size r , and 
propagating as fronts of light waves. Regions (a) and (b) represent their internal regions, 
where the stress tensor has identical forms, (|2.7|) . In the external regions (c), (d) the 



metric ( |2.3|) changes to 

ds 2 = (r 2 /r 2 - \y l dr 2 - (r 2 /r 2 - l)dt 2 - r 2 (d6 2 + sin 2 9d<p 2 ) , (2.10) 

and describes an expanding (or shrinking) universe where r now plays the role of time and 
t has the meaning of a spatial coordinate. In these case, the horizon surfaces are removed 
to spatial infinities (t = ±00), and the structure of the stress tensor can be found on 



the base of (|2.7|) . It follows immediately that in these domains T"((3) corresponds to a 
non equilibrium thermal system. Although the energy flows from the infinities (t = ±00) 
are absent, the energy density decreases to the vacuum one due to the factor [<?tt(r)]~ 2 in 
according with adiabatic expansion of this universe. 

A remark concerning the thermodynamical parameters of the bubbles like their energy, 
and entropy it is worth being done here. These quantities can be computed in two ways: 
through the local energy Tf(/3) [fj, or in terms of the one-loop partition function Q. The 
first method cannot be straightforwardly applied due to the divergence of the local energy 
and entropy at the horizon, whereas the other one gives a finite answer under (^-function 
regularization [|J . This fact might be understood in terms of renormalization of bubbles 
surface energy, which can cancel the local divergence. To partially support the above 



statement we refer to the analysis carried out for manifolds with conical singularities [ 14] . 
Such kind of defects is inherent to "gravitational instantons" associated to static de Sitter 
space (U) at p H M- 
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3 The semiclassical analysis 

Consider now changing of de Sitter geometry induced by the vacuum polarization in the 
given quantum states. Einstein equations in semiclassical approximation are known to 
read as 

R^{g) - ^R(g)g»u + + a ^H^g) + (3 {2) H^(g) = -8nG < % v > p . (3.11) 



Here g^ u = g^ v + h^ u indicates the solution of (|3.11| ) including quantum corrections h^ u 
to de Sitter metric g^, and A = 3/r 2 is the cosmological constant. As for the additional 
tensors ^'H^ and ^H^ u that appear due to the quantum corrections, and are defined 
in terms of the geometrical quantities ||, we will neglect them in the further analysis 
assuming the unknown constants a, (3 in ( |3.11| ) to be equal to zero. Then, expanding 



3.11 ) up to the first order in the metric perturbation h^ v , we have 



r \ r / 

(3.12) 

An essential progress in solving these equations can be achieved observing that the 
anomalous part of the renormalized stress tensor in the r.h.s. of ( [3.12|) only gives rise to 
a redefinition of the de Sitter radius tq. Taking into account this trivial effect we can use 
instead of the total tensor only its traceless part T"{/3). Thus, after a little transformation, 
the initial equations take the form 

^ + lV = - 16ffG W) > ( 3 - 13 ) 
r o 

where we imposed the gauge conditions h j v = 0, M = (the latter is compatible with 
the dynamical equations and fixes the residual gauge freedom). Besides this, the property 
of de Sitter space for which R\ vm = TQ 2 {gx ll g va — g\ a gvn) has been used in ( |3.13| ). 

To proceed in our analysis, it is worth reminding that in the given quantum state 
(P 7^ Ph) the total de Sitter symmetry 50(1,4) is broken to the subgroup T\ x 50(3), 
where T x stands for de Sitter boosts associated with translations along the coordinate t. 
This enables us to restrict the degrees of freedom of h^, to the only nonzero components 
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htt(r), h rr (r), h tr {r) = h rt (r), and h ik = — a{r)r 2 ^ ikl where the indexes i, k are referred to 



the 9, ip coordinates, and 7^ is the metric on S 2 (in accordance with 



9%k = -r'~fik)- 



Due to the residual symmetry and gauge conditions imposed, there is only one indepen- 



dent equation among ( |3.13|) . After some algebra one can show that in static coordinates 
it reduces to a second order differential equation for h rr 



4(1-Z/) J 2 



d 2 _ 26y - 16 d _ 28y - 8 
dy 2 y dy y 2 



h rr (y) 



(3.14) 



written in terms of the variable y(r) = 1 — r /r 2 and a constant 



K((3) 



90^ 



(3.15) 



Remarkably, the homogeneous equation associated to ( J3.14 ) admits two simple solutions: 
y~ 2 , and [y 2 (l — y) 3 ^ 2 ]^ 1 . Thus the general integral of ([3.14D can be given in the following 
form 



h rr (r) 



+ 



4y 2 (r) 
A 



2 2 1 

3 + l-y(r) + [l-y(r)f/ 2 

B 



log 



l + Jl-y(r) 



(3.16) 



y 2 { r ) ?/ 2 ( r )[l — y(r)} 3 / 2 
It is worth mentioning that coordinate systems in de Sitter spaces, (|2.3|), fl2.10| ), can be 
treated on equal footing if one adds a small imaginary part to de Sitter radius r — > r + ie. 
This regularization preserves the structure of all the equations in both quantum and 
classical theory. It removes the singularity at the horizons in such a way that ( |2.1U|) , and 
( [2.3|) can be unified in a single expression for the metric valid for < r < 00. However, 
although such regularization allows to define the solution in all the region < r < 00, the 
integration constants A and B cannot be chosen to be equal in all this interval. Indeed, 
if A and B are fixed in some way in the inner domain (a), then after passing in the 
external region, say (c), one obtains due to the logarithm complex values for the function 
h rr (r). It means that integration constants for external and internal problems in the given 
semiclassical approach should be found independently, on the base of additional physical 
motivations. 

In the static regions A and B can be chosen from the condition h^ u {r) = at r = that 
corresponds to the natural assumption for the vacuum effects to disappear in the flat limit 



r/r << 1, when the background space converts into Minkowski one. As for the external 
domains, one can require for the perturbed metric to approach the de Sitter metric in 
the limit r — > oo, where local quantum effects disappear (see fl2.7|) ). The components h tt , 
hik, can be obtained from (|3.16|) using constraints. On the other hand, the non diagonal 
element h tr , fixed by the condition h tu ,u = 0, turns out to be completely independent of 
the stress tensor and other components. For this reason, we put it to be zero everywhere, 
which inside the horizon can be also justified by claiming its regularity at r = 0. 

By virtue of the expression ( |3.16| ) and conditions chosen we can represent all the non 
zero components of the metric perturbations in the following form 



h tt {r) 



K{fl) 



4r 2 r fr -r 

2 -^ + ^ log UT7 



+ A 3 



2r 2 



+ B- 



r 



(3.17) 



h rr (r) 



K{P) 
4 y 2 (r) 



2r 



r 



r + r 



A 



B- 



y 2 (r) y 2 (r)r 3 



(3.18) 



a{r) 



K(P) 



o 



+ ^log 



_2 + ^ 
3 r 2 y(r) ' 4r 3 



r — r 



2r 3 



(3.19) 



>r + r, 

Following the previously mentioned requirements, the integration constants A and 
B have to be put equal to zero for < r < r , whereas for r < r < oo we have 
A = —K{/3)/6. As far as the value of B in the external region is concerned, it is left 
undetermined. This fact can be understood reminding that in the approximated approach 
of semiclassical theory two regions are separated by an infinite barrier, whereas in the 
correct theory both solutions should be matched at the horizon, so fixing the remaining 
constant. Note also that solutions ( |3.17|) -( |3.19| ) are valid only in the region where the 
perturbation expansion is reliable. In terms of the Planck length Ipi = vG, this region 
results to be \r — r | >> |1 — /3^//3 4 | 1//2 /p/. Therefore, if the de Sitter radius r ~ Ipi the 
present analysis is not applicable to the bubble interior. 

Finally, we observe that solutions ( |3.17 )-( 3.19p of semiclassical Einstein equations in de 
Sitter background, can be immediately extended to anti-de Sitter space just substituting 
ro with itq. This is connected to the fact that the renormalized stress tensor at finite 
temperature in anti-de Sitter space has the same structure of de Sitter one, once a complex 
scalar field with Neumann and Dirichlet boundary conditions for its two independent 
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components is chosen 







4 Gravitational effects 



We discuss now the properties of the space taking into account quantum corrections ( 3.17|) 
- ( p,19| ) to de Sitter metric. It can be done considering the different models of de Sitter 
space corresponding to flat, open or closed expanding universes. For sake of simplicity we 
restrict in this letter our analysis to the flat model, which is generally used for cosmological 



applications. A more complete study will be given in a forthcoming paper [ IB| . 
The corresponding metric has the familiar form 

ds 2 = dr 2 - e 2r/ro (d£ 2 + f W) , < £ < oo , (4.20) 

(dQ 2 = d6 2 + sin 2 ^^ 2 is the line element on a sphere). This system of coordinates 
maps only a half of the de Sitter hyperboloid. For a more clear representation we assume 
it to cover regions (a) and (c) on figure 1. In this case, one of the singularities of the 
global stress tensor is placed inside the universe (}4.20|), in such a way that bubble's center 



corresponds to the origin of coordinates. At the same time, the other one is at spatial 
infinity (£ — > oo). The map between Q4.20|) and ( |2.3| ) or (|2.I0|) can be written as follows 



t = r-^log(l-^e 2 ^J , r = £e r/ro . ( 42 r 



Note that by definition, the shell of the bubble coincides with de Sitter horizon thus its 
size remains unchanged during inflation 3 , but due to the expansion every point, once 
inside it, will move away. Near the matter distribution surface £o = r exp(— r/r ) the 
metric (|4.20|) changes to 

ds 2 = (1 + h m ) drf + 2h^ drjd^ - [e 2ri/r ° - d£ 2 - e 2r?/r °£ 2 (1 + a) dVt 2 , (4.22) 

where we introduced rj instead of r to indicate that now it is not the proper time. The 
metric perturbations h vv , h v ^ and follow from ( |3.17| ), ( |3.18| ) and map ( f4.21| ), where r 
is replaced by r/ 

M^,0 = - 2 h tt + %e 2 ^h rr , (4.23) 

y n 



3 A discussion about this point is given in section 5. 



1 



r 



h tt + K 



3 2??/ro 



(^L e lr,/ro htt + h \ 

W y J 



(4.24) 



(4.25) 



In Eqs. ( |122j )-( |4T2"5| ) all the functions y, a, h a and h rr are supposed to be expressed in 
terms of rj and £. So far as rj does not represent the proper time we need an additional 
transformation to pass to comoving coordinates (r,x,9,(p), similar to ( |4.20| ), 



ds > = dT 2 _ e 2r/r 



(1 + h xx (x,T))dx 2 + x 2 (l + h nn (x,r))dQ 2 = dr 2 - -f ij dx i dx j (4.26) 



with i,j = 1,2,3, corresponding to x, 9 and ip respectively. The above transformation 
reads 



V = r + r]i(r,x) 



(4.27) 



and the functions rji and £i obey to the following differential equations 



2^ + h m = 



e 2 ^-h r dVl 







dr ^ dx 

In terms of these functions the metric corrections h xx and result to be 



(4.28) 



l>r. ■ =2^ + 2^i-e- 2 ^% , h nn = 2^ 
r ox r 



h 



£i , 
2 h a 

x 



(4.29) 



It is worth observing that the integration constant for £i, obtained from ( |4.28| ), is not 
relevant because it can always be changed by redefinition of x coordinate in ( 4.26Q . As far 
as 771 is concerned, its value is fixed, up to an inessential numerical constant, by requiring 
that metric components ( |4.26| ) depend only on the distance r = xexp(r/ro) from the 
center of the perturbation. 



The metric components can be readily obtained from ( |4.28|) on the total space (0 < 
r < r and r < r < oo). For sake of brevity we only present the expressions for the 
outside region to investigate the metric asymptotic at infinity. They read 



h 



K{fi) 



log 



4 ° \ r 2 — Tn 



+ 7T^ log 



r — r 
r + r 



+ 



r 2 2r 2 - 3r 2 



(4.30) 
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h xx = K(P) 



5 , ( r 2 \ Tq r 2 — 3r 5 



-l^gi——, + 



4 y r 2 _ r 2 J ^ r 2 r 2 _ r 2 



(4.31) 



where the undefined constant B was neglected, because it does not affect the long-range 
behavior (r ^> r ) of perturbations. By using ( [4.30| ) and ( |4.31| ) one can see that gravita- 
tional field far from the bubbles decreases like r~ 2 . It corresponds, due to the expansion, 
to an exponential falling with time r for a point of given x. 

The effects of the quantum corrections to the metric ( |4.20| ) can be expressed in another 
way by computing the spatial curvature (7) of the universe Q4.26|) 

^ (3) (7) = /'•;! - h$ 1 = |e- 2 ^ ^(x/(x, r)) , (4.32) 

d d 
f(x, t) = —{x haa) - h xx = 3a + x—a . (4.33) 
ax ax 

Remarkably, due to (f03|) it turns out to be completely independent of the constant B 
and takes the simple analytical forms 

Qiv«2 ^,2 ^,2 ^2 I ^-.2 

tf) WlK r » = ^)^ , R {3) ^)\r>r =K(P)f 2V ^ 2 , (4.34) 



for inside and outside regions respectively. By definition ( |3.15| ) of K(j3) it follows that 



curvature is negative if the bubble temperature [3~ l higher then the Hawking one Z^ 1 , 
and changes the sign when < fijj . Moreover, its value at r = does not depend 
on time i?( 3 )(7)| r=0 = 3K(fi)rQ 2 , whereas in external region decreases as r -4 . This 
long range behaviour of the gravitational field is connected to the massless nature of the 
chosen scalar field. For massive matter fields one can expect that the gravitational effects 
will be exponentially smeared, on a length of the order of the inverse mass, outside the 
bubble shell. 

5 Conclusions and remarks 

In this letter the local properties of thermal states for a quantum field in de Sitter space- 
time have been investigated. We interpret the known singularity of renormalized stress- 
tensor as a matter distribution located on two spherical surfaces (boundaries of bubbles) 
moving far away one from each other. Such domains, characterized by temperature, 
entropy, and etc., give rise to a number of semiclassical effects. 
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To find the gravitational field produced by these perturbations, the linearized Einstein 
equations on de Sitter background have been exactly solved for stress tensor of a conformal 
scalar field at finite temperature. The computations show the gravitational effects are 
weak far from the bubbles, and in the expanding universe they are exponentially damped. 

In realistic cosmological models both the radius and thermodynamical properties of 
the bubbles might change during the universe evolution. Note that external region of 
static coordinates ( |2.3|) can be considered as a "black hole" absorbing "information" 
from the inner domain due to decays there of unstable field configurations of the Higgs 
scalars. As a result of this process, the bubble size, associated to the horizon, grows 
up and achieves cosmological values determined by the inverse of the Hubble constant, 
H^ 1 m 10 28 cm at the present. However, due to the exponential damping, the presence of 
these inhomogeneous regions does not seem to affect the observable part of the universe, 
provided if it is far enough from them. We confine the present letter to these short 
remarks, leaving more complete analysis of the cosmological aspects of such theory to a 
forthcoming paper [|l(| . 

Finally, it is worth observing that the solutions obtained for the semiclassical Einstein 
equations in de Sitter background can be straightforwardly extended to the anti-de Sitter 
geometry. In this case in fact, the stress tensor has the same structure of the de Sitter 
one, once a complex scalar field with appropriate boundary conditions is chosen |L5] . 
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Figure caption 

Figure 1: In this diagram left and right edges have to be identified. Regions (a) and 
(b) correspond to static coordinates with < r < r , separated by light cones from (c) 
and (d) regions where r < r < +00. Dashed lines in regions (a) and (b) represent the 
trajectories of a particle and its antipod being in the coordinate origins. 
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